Consistent Testing for Recurrent Genomic Aberrations

Genomic aberrations, such as somatic copy number alterations, are frequently observed in tumor tissue. Recurrent aberrations, occurring in the same region across multiple subjects, are of interest because they may highlight genes associated with tumor development or progression. A number of tools have been proposed to assess the statistical significance of recurrent DNA copy number aberrations, but their statistical properties have not been carefully studied. Cyclic shift testing, a permutation procedure using independent random shifts of genomic marker observations on the genome, has been proposed to identify recurrent aberrations, and is potentially useful for a wider variety of purposes, including identifying regions with methylation aberrations or overrepresented in disease association studies. For data following a countable-state Markov model, we prove the asymptotic validity of cyclic shift $p$-values under a fixed sample size regime as the number of observed markers tends to infinity. We illustrate cyclic shift testing for a variety of data types, producing biologically relevant findings for three publicly available datasets.Comment: 35 pages, 7 figure

DiNAMIC - A Method for Assessing the Statistical Signficance of DNA Copy Number Aberrations

DNA copy number gains and losses are commonly found in tumor tissue, and some of these aberrations play a role in tumor genesis and development. Although high resolution DNA copy number data can be obtained using array-based techniques, no single method is widely used to distinguish between recurrent and sporadic copy number aberrations. Here we introduce Discovering Copy Number Aberrations Manifested In Cancer (DiNAMIC), a novel method for assessing the statistical signifcance of recurrent copy number aberrations. DiNAMIC uses two resampling schemes - a permutation method and a bootstrap procedure - both of which largely preserve the correlation structure found in the underlying DNA copy number data. It is important to maintain as much of the correlation structure as possible when resampling, and we believe this may yield additional power to detect recurrent aberrations. Extensive simulation studies show that DiNAMIC controls false positive discoveries in a variety of realistic scenarios. We use DiNAMIC to analyze two publicly available tumor datasets, and our results show that DiNAMIC detects multiple loci that have biological relevance. Although DiNAMIC provides methods for detecting CNAs, loci that exhibit aberrant copy number do not always lie in genes related to the tumor phenotype. Therefore we introduce methods for computing confidence intervals around CNAs. Copy number datasets often contain data obtained from subjects with different subtypes of a given tumor type, and the tumor subtypes can harbor distinct genetic mutations. Because groups of subjects may have similar copy number probles, we present methods for determining which subjects contribute to a given CNA. Some studies collect both DNA copy number and clinical data from subjects, but no methods for jointly analyzing both data types are widely used. We describe a preliminary testing procedure for comparing the locations of CNAs and loci whose copy number values are highly associated with a given clinical variable

Hedgehog-GLI Signaling Inhibition Suppresses Tumor Growth in Squamous Lung Cancer

Lung squamous cell carcinoma (LSCC) currently lacks effective targeted therapies. Previous studies reported overexpression of HEDGEHOG (HH)-GLI signaling components in LSCC. However, they addressed neither the tumor heterogeneity nor the requirement for HH-GLI signaling. Here, we investigated the role of HH-GLI signaling in LSCC, and studied the therapeutic potential of HH-GLI suppression

Homozygosity for the SCN10A Polymorphism rs6795970 Is Associated With Hypoalgesic Inflammatory Bowel Disease Phenotype

Background: Hypoalgesic inflammatory bowel disease (IBD), a condition in which patients with active disease do not perceive and/or report abdominal pain, is associated with serious complications and there is a lack of cost-effective, reliable diagnostic methods to identify “at-risk” patients. The voltage-gated sodium channels (VGSC's), Nav1.7, Nav1.8, and Nav1.9, are preferentially expressed on nociceptive neurons, and have been implicated in visceral inflammatory pain. At least 29 VGSC single nucleotide polymorphisms (SNPs) have been implicated in chronic somatic pain syndromes, but little is known about their role in human visceral sensation. We hypothesized that disruptive VGSC polymorphisms result in anti-nociceptive behavior in IBD.Methods and Findings: We performed targeted exome sequencing and/or TaqMan genotyping to evaluate the Nav1.7, Nav1.8, and Nav1.9 genes (SCN9A, SCN10A and SCN11A) in 121 IBD patients (including 41 “hypoalgesic” IBD patients) and 86 healthy controls. Allelic and genotypic frequencies of polymorphisms were compared among study groups who had undergone characterization of intestinal inflammatory status and abdominal pain experience. Forty-nine total exonic SNPs were identified. The allelic frequency of only one non-synonymous SNP (rs6795970 [SCN10A]) approached significance in hypoalgesic IBD patients when compared to other IBD patients (p = 0.096, Fisher's exact test). Hypoalgesic IBD patients were more likely to be homozygous for this polymorphism (46 vs. 22%, p = 0.01, Fisher's exact test).Conclusions: This is the first human study to demonstrate a link between a genetic variant of SCN10A and abdominal pain perception in IBD. These findings provide key insights into visceral nociceptive physiology and new diagnostic and therapeutic targets to consider in IBD and other gastrointestinal conditions associated with chronic abdominal pain. Further studies are required to elucidate the precise pathophysiological impact of the rs6795970 polymorphism on human gastrointestinal nociception

Integrated RNA and DNA sequencing improves mutation detection in low purity tumors

Identifying somatic mutations is critical for cancer genome characterization and for prioritizing patient treatment. DNA whole exome sequencing (DNA-WES) is currently the most popular technology; however, this yields low sensitivity in low purity tumors. RNA sequencing (RNA-seq) covers the expressed exome with depth proportional to expression. We hypothesized that integrating DNA-WES and RNA-seq would enable superior mutation detection versus DNA-WES alone. We developed a first-of-its-kind method, called UNCeqR, that detects somatic mutations by integrating patient-matched RNA-seq and DNA-WES. In simulation, the integrated DNA and RNA model outperformed the DNA-WES only model. Validation by patient-matched whole genome sequencing demonstrated superior performance of the integrated model over DNA-WES only models, including a published method and published mutation profiles. Genome-wide mutational analysis of breast and lung cancer cohorts (n = 871) revealed remarkable tumor genomics properties. Low purity tumors experienced the largest gains in mutation detection by integrating RNA-seq and DNA-WES. RNA provided greater mutation signal than DNA in expressed mutations. Compared to earlier studies on this cohort, UNCeqR increased mutation rates of driver and therapeutically targeted genes (e.g. PIK3CA, ERBB2 and FGFR2). In summary, integrating RNA-seq with DNA-WES increases mutation detection performance, especially for low purity tumors

<html>Autologous reconstitution of human cancer and immune system <i>in vivo</i></html>

Correlative studies from checkpoint inhibitor trials have indicated that better understanding of human leukocytic trafficking into the human tumor microenvironment can expedite the translation of future immune-oncologic agents. In order to directly characterize signaling pathways that can regulate human leukocytic trafficking into the tumor, we have developed a completely autologous xenotransplantation method to reconstitute the human tumor immune microenvironment in vivo. We were able to genetically mark the engrafted CD34+ bone marrow cells as well as the tumor cells, and follow the endogenous leukocytic infiltration into the autologous tumor. To investigate human tumor intrinsic factors that can potentially regulate the immune cells in our system, we silenced STAT3 signaling in the tumor compartment. As expected, STAT3 signaling suppression in the tumor compartment in these autologously reconstituted humanized mice showed increased tumor infiltrating lymphocytes and reduction of arginase-1 in the stroma, which were associated with slower tumor growth rate. We also used this novel system to characterize human myeloid suppressor cells as well as to screen novel agents that can alter endogenous leukocytic infiltration into the tumor. Taken together, we present a valuable method to study individualized human tumor microenvironments in vivo without confounding allogeneic responses

Genetic evolution of uveal melanoma guides the development of an inflammatory microenvironment

Experimental cancer immunology and therap

Cancer-Derived Mutations in KEAP1 Impair NRF2 Degradation but not Ubiquitination

NRF2 is a transcription factor that mediates stress responses. Oncogenic mutations in NRF2 localize to one of its two binding interfaces with KEAP1, an E3 ubiquitin ligase that promotes proteasome-dependent degradation of NRF2. Somatic mutations in KEAP1 occur commonly in human cancer, where KEAP1 may function as a tumor suppressor. These mutations distribute throughout the KEAP1 protein but little is known about their functional impact. In this study, we characterized 18 KEAP1 mutations defined in a lung squamous cell carcinoma tumor set. Four mutations behaved as wild-type KEAP1, thus are likely passenger events. R554Q, W544C, N469fs, P318fs, and G333C mutations attenuated binding and suppression of NRF2 activity. The remaining mutations exhibited hypomorphic suppression of NRF2, binding both NRF2 and CUL3. Proteomic analysis revealed that the R320Q, R470C, G423V, D422N, G186R, S243C, and V155F mutations augmented the binding of KEAP1 and NRF2. Intriguingly, these 'super-binder' mutants exhibited reduced degradation of NRF2. Cell-based and in vitro biochemical analyses demonstrated that despite its inability to suppress NRF2 activity, the R320Q 'superbinder' mutant maintained the ability to ubiquitinate NRF2. These data strengthen the genetic interactions between KEAP1 and NRF2 in cancer and provide new insight into KEAP1 mechanics

BRG1/SMARCA4 Inactivation Promotes Non-Small Cell Lung Cancer Aggressiveness by Altering Chromatin Organization

SWI/SNF chromatin remodeling complexes regulate critical cellular processes including cell cycle control, programmed cell death, differentiation, genomic instability and DNA repair. Inactivation of this class of chromatin remodeling complex has been associated with a variety of malignancies, including lung, ovarian, renal, liver and pediatric cancers. In particular, ~10% of primary human lung non-small lung cancers (NSCLC) display attenuations in the BRG1 ATPase, a core factor in SWI/SNF complexes. To evaluate the role of BRG1 attenuation in NSCLC development, we examined the effect of BRG1 silencing in primary and established human NSCLC cells. BRG1 loss altered cellular morphology and increased tumorigenic potential. Gene expression analyses showed reduced expression of genes known to be associated with progression of human NSCLC. We demonstrated that BRG1 losses in NSCLC cells were associated with variations in chromatin structure, including differences in nucleosome positioning and occupancy surrounding transcriptional start sites of disease-relevant genes. Our results offer direct evidence that BRG1 attenuation contributes to NSCLC aggressiveness by altering nucleosome positioning at a wide range of genes, including key cancer-associated genes